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We have studied the magnetization of various, well characterized samples of highly oriented 
pyrolitic graphite (HOPG), Kish graphite and natural graphite to investigate the recently reported 
ferromagnetic-like signal and its possible relation to ferromagnetic impurities. The magnetization 
results obtained for HOPG samples for applied fields parallel to the graphene layers - to minimize 
the diamagnetic background - show no correlation with the magnetic impurity concentration. Our 
overall results suggest an intrinsic origin for the ferromagnetism found in graphite. We discuss 
possible origins of the ferromagnetic signal. 



INTRODUCTION 

Wc find in the literature reports on ferromagnetism in 
carbon-based compounds which apparently did not at- 
tract the necessary attention in the community. For ex- 
ample there are reports on ferromagnetic-like signals in 
pyrolytic carbon j[| and in HOPG [g. In the first study 
pyrolytic carbon with higher crystallinity than usual was 
prepared from adamantane by chemical vapor deposition 
0. The metal impurities were claimed to be < 25 ppm 
and the ferromagnetic-like loop was reported to persist 
up to 400 K. In other work || the main study was 
the magnetization of Pd nanoparticles encapsulated in a 
graphite host. In the inset of their Fig. 2 the authors 
show a ferromagnetic-like hysteresis loop for graphite 
without Pd nanoparticles at a temperature of 300 K. 
We note that the saturation magnetization M s ~ 10~ 2 
emu/g is much larger than the one measured in this work 
and by Kopelevich et al. j|. Because the impurity con- 
centration of that sample was not reported the origin 
of this signal remains, however, unclear. Amorphous- 
like carbon prepared by direct pyrolysis was found to be 
a stable strong magnet with a saturation magnetization 
comparable to magnetite The effect of the pyrolysis 
temperature on similar materials was also reported [^J. 
Ferromagnetism in pure C-60 induced by photo-assisted 
oxidation was reported to remain up to a critical tem- 
perature of 800 K, although the authors stated that only 



a small part of the sample contributes to the ferromag- 
netism [H . More recently, ferromagnetism in pure rhom- 
bohedral polymerized C-60 samples was reported to re- 
main up to 500 K @. 

The recent report of ferromagnetic- and 
superconducting-like H magnetization hysteresis loops 
in HOPG in a broad temperature range opens several 
questions about the origin of the observed anomalies. 
Kopelevich et al. |}| found that the magnetization of 
HOPG samples shows ferromagnetic-like hysteresis loops 
up to 800 K. The details of the hysteresis depend on the 
sample, sample heat treatment and on the direction of 
the applied field. There has been a controversy between 
the opinion of the referees on the possible influence 
of impurities on the observations Therefore, more 
systematic work on the observed signal was necessary. 

As clearly stated by J.S. Miller |s) for magnetic 
polymers, it appears that most of the claimed high- 
temperature ferromagnetism in nominally magnetic-ion 
free compounds in the past turned out to be of extrinsic 
origin. To establish that the low saturation magnetiza- 
tion of a material, as in the case of the reported HOPG 
(M s < 3mG = 3 x 10 -3 emu/cm 3 ) might be of intrinsic 
origin, reproducibility and a quantification of extrinsic 
contributions must be achieved. Therefore and in order 
to clarify whether the observed ferromagnetism found in 
Ref. || reflects an intrinsic property of graphite or is due 
to ferromagnetic impurities, we have measured different 
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oriented graphite samples with different magnetic impu- 
rity concentrations (mainly Fe). 

Both ferromagnetism and superconductivity have been 
predicted to arise in graphene layers in regions of certain 
topological defect structures j| . On the other hand, fer- 
rimagnetic behavior due to edge states is also expected 
in graphene layers j[o|. Taking into account all the ex- 
perimental and theoretical work cited above we feel that 
a systematic study of the ferromagnetic-like hysteresis 
loops found in HOPG || as a function of magnetic impu- 
rities is necessary before speculating on its origin. Only 
a careful analysis of the data accompanied by a complete 
impurity concentration study on the same samples may 
provide an answer on the origin of the ferromagnetism 
found in HOPG. 

The reported ferromagnetic-like signal was measured 
mostly for an applied field in direction parallel to the 
planes, whereas superconducting-like loops in HOPG 
were mainly observed with the magnetic field applied 
perpendicular to the planes after subtraction of the dia- 
magnetic background || . This diamagnetic contribution 
depends strongly on the angle between applied field and 
the graphene layers reaching a minimum for the paral- 
lel orientation. Therefore, and in order to minimize any 
possible error due to the subtraction of the diamagnetic 
contribution, all the results reported in the current study 
were obtained with the field parallel to the graphene lay- 
ers. With respect to the superconducting-like behav- 
ior observed in graphite || we would like to note that 
clear evidence for superconductivity in simple mixtures 
of graphite powder with sulphur has been reported re- 
cently 0. 

The paper is organized as follows. In the next Section 
the samples characteristics and experimental details are 
given. In Section 3 we describe the results and in Section 
4 we discuss the observed behavior and its possible origin. 
A short summary is given in Section 5. 

SAMPLE CHARACTERISTICS AND 
EXPERIMENTAL DETAILS 

Table I shows the magnetization characteristics, the 
impurity concentrations, the masses as well as the full 
width at half maximum (FWHM) of x-ray rocking curves 
in 9 (AO) and in 29 (A(29)) of the measured graphite 
samples. Samples HOPG-2 and -3 were prepared at the 
Research Institute "Graphite" (Moscow) as described in 
Ref. Two commercial HOPG samples were from 

Union Carbide (samples UC3 and UC4, cut from the 
same batch) and Advanced Ceramics (samples AC1, 
batch 10158, and AC2, batch 25672). Furthermore, we 
have characterized a Kish and a natural graphite crystal, 
this last obtained from Ticonderoga (USA) |T^| . 

The magnetization measurements were done with a 
Superconducting Quantum Interference Device (SQUID) 



from Quantum Design with a 7 T superconducting 
solenoid. In general we have selected three cm scan 
length for the magnetization measurements. The SQUID 
response enables measurements of samples with a mini- 
mum magnetic moment m~3x 1CU 10 Am 2 (3 x 10~ 7 
emu) . 

Because the magnetization signals of the samples are 
rather small we have checked the background signal from 
the sample holders in the SQUID. This contribution 
turned out to be non-negligible, and special care has been 
applied in the selection of the holder as well as with the 
sample handling procedure in order to prevent contami- 
nation of the holders and/or samples. For measurements 
below 350 K we have used a plastic tube with a small 
window and an adhesive tape (Tesa AG, Hamburg). The 
background signal of this sample holder with tape but 
without sample is shown in Fig. 1. Small but measurable 
hysteresis is observed which lies close to the resolution 
limit of our SQUID magnetometer. For measurements 
above 350 K the samples were fixed in a quartz tube be- 
tween two quartz spacers and enclosed under a pressure 
of 2 x 10 -6 torr. The background signal in this case was 
negligible in comparison with the sample response. 

In general, the ferromagnetic hysteresis is superim- 
posed to a linear in field diamagnetic background. As 
noted above, this diamagnetic signal depends on the mis- 
alignment of the sample with respect to the field as well 
as on the intrinsic misalignment of the graphene planes in 
the sample. The diamagnetic contribution for each sam- 
ple has been determined in-situ and in a straightforward 
way: Using a linear regression through the data points 
obtained for magnetic fields up to 5 T we obtained the 
diamagnetic slope x\\ = M(H)/H. We note that the 
magnetization due to the superimposed ferromagnetic 
signal saturates at fields larger than 2 kOe and there- 
fore it does not influence the linear field dependence of 
the diamagnetic background. The obtained susceptibil- 
ity ranges between > X||(10i^) > —5.3 x 10~ 6 emu/gOe 
for the studied samples. The ferromagnetic loop is then 
obtained from the measured magnetization M(H) sub- 
tracting the diamagnetic background as M(H) —x\\H. 

The crystallinity of the graphite samples was investi- 
gated by high resolution x-ray diffraction with a Philips 
X'pert diffractometer. 9 and 29 scans around the sym- 
metrical (002) reflection were carried out using the Cu- 
Kal radiation. Furthermore, for three samples we have 
determined the average in-plane correlation length D 
from the 29 scans around the symmetrical (110) reflec- 
tion. As expected, D correlates with the FWHM ob- 
tained from the 9 (002)-scans. Table I shows the FWHM 
of the 110 reflection as well as D for the samples AC2, 
HOPG-2 and HOPG-3. 

The measurements of the impurity concentration were 
carried out with Particle Induced X-ray Emission (PIXE) 
using a 2 MeV proton beam. For the samples studied 
here the typical minimum detection limit (MDL) of this 
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method was < 2/ig/g (for some elements, much smaller). 
With exception of the natural graphite sample, the im- 
purity concentration was distributed homogeneously in 
the samples. The concentration of all other magnetic 
elements not listed in Table I was less than 1 /xg/g. 

Scanning tunelling microscopy (STM) was performed 
with a DI Nanoscope E in conjunction with a Small Sam- 
ple scan head ("A" scanner, Digital Instruments, Santa 
Barbara, CA). Pt/Ir (80:20) tunneling tips (wire diam- 
eter: 0.25 mm) were from Piano (Wetzlar, Germany). 
Images were obtained under ambient conditions with a 
tip bias of typically +20 mV in the constant current- 
mode (scan rate / = 3 ... 6 Hz; current Iq rj 2 nA) or 
in the constant height-mode (scan rate / = 20 ... 60 Hz). 
While in both modes, a clear-cut discrimination of topol- 
ogy (height) differences from density-of-state (current) 
differences cannot be afforded, the first scanning mode 
emphasizes the topological structure whereas the latter 
accentuates local variations of the density of states. 

RESULTS 

Figures |^ to |] show examples of the measured mag- 
netization loops obtained in some of the samples at 300 
K and 10 K. All HOPG samples show ferromagnetic-like 
hysteresis of the same order of magnitude. Figure ^] shows 
the hysteresis loops obtained for the samples from Ad- 
vanced Ceramics at two different temperatures. Whereas 
the loop shown in (b) is obtained after subtraction of a 
diamagnetic background signal X|| = — 1.32 x 10~ 6 emu/g 
Oe, the loop in (a) was obtained directly from the mea- 
surements without subtraction of any diamagnetic back- 
ground. This result and the range of values obtained for 
X\\ suggest that the diamagnetic signal measured at high 
fields applied nominally parallel to the planes of HOPG 
samples is strongly influenced by the misalignment of the 
sample. 

Figure |^ shows the loops obtained for two Union Car- 
bide samples of different masses cut from the same batch. 
These results provide a measure of the typical scatter of 
the main values of the hysteresis loops in similar samples. 

In order to check if there is any correlation of the 
hysteresis values with the measured Fe-concentration we 
plot in Figs. [|, ^ and [?] the remanent magnetization M ri 
the magnetization at 2kOe M(2kOe) and the coercive 
field H Cl respectively, as a function of impurity content. 
From Fig. |5| one would tend to infer a weak increase 
of M r with Fe concentration. However, the remanent 
magnetization is not the characteristic quantity that one 
expects to change systematically with Fe concentration 
since it depends primarily on the pinning of magnetic 
domains. Rather, a correlation of the saturation magne- 
tization with the impurity content would be significant. 
Figure || shows that such a correlation is not observed 
within experimental error. A similar result is obtained 



for the coercive field, see Fig. |7]. 

The FWHM obtained in 8 and 29 scans provides a mea- 
sure of the variation in the orientation of a particular set 
of lattice planes (mosaic spread) and in the spacing of the 
graphene planes (e.g., inferred by dislocations), respec- 
tively. Figure || shows the magnetization at saturation as 
a function of A8 (a) and A (29) (b). We find no system- 
atic dependence of M(2kOe) on A9. A clear increase of 
M(2kOe) with A(29) can be recognized in Fig. |(b). 

It has been reported H that annealing of HOPG sam- 
ples leads to the disappearance of the superconducting- 
like hysteresis loop and turns it into a ferromagnetic-like 
loop. Also, short time annealing (~ 2 h) increases the 
ferromagnetic-like signal . In the current work we have 
measured the hysteresis as a function of the annealing 
time and temperature. An increase of the ferromagnetic 
signal at the beginning of the annealing procedure was 
also observed in this work. Longer time annealing re- 
duces however the ferromagnetic signal and reaches sim- 
ilar values as before the annealing. In Fig. [|(b) we show 
the data obtained at 300 K after annealing a different 
piece of the HOPG-2 sample for 16 h at 700 K (o) and 
additionally after 19 h at 800 K (a). The annealing pro- 
cedure was done in situ in the SQUID under vacuum. It 
is observed that long time annealing has a weak influ- 
ence, if any, on the hysteresis width and the saturation 
magnetization. 

DISCUSSION 

The origin of the ferromagnetic-like signal in HOPG 
samples is controversial. Therefore, we will discuss below 
different possible origins in view of the existing evidence. 

(a) Magnetic impurities. The possibility of a signifi- 
cant contribution of magnetic impurities to the measured 
magnetic properties is taken serious and carefully ana- 
lyzed before we discuss possible intrinsic origins of the 
hysteresis loops measured in graphite. In spite of the fact 
that no correlation between the observed magnetization 
and the Fe-concentration exists, one might still believe 
that those impurities are sufficient to produce the mea- 
sured signals. A typical argument used to attribute a 
ferromagnetic signal to impurities would be as follows: 
Let us assume that all Fe impurities would form iron or 
magnetite (Fe3 04) clusters and that all clusters are large 
enough to behave ferro- or ferrimagnetically. An impu- 
rity content of 1 fj,g Fe per gram graphite would then 
contribute ~ 2.2 x 10~ 4 emu/g to the magnetization in 
the case of Fe clusters or ~ 1.4 x 10~ 4 emu/g in the case 
of Fes04 clusters assuming that these clusters behave 
like bulk materials. Considering a linear increase of the 
magnetic moment with the Fe concentration, for most of 
the samples the possible contribution of the Fe impuri- 
ties exceeds clearly the measured magnetization values 
at saturation or are on the same order of magnitude, see 
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Fig. [| In the case of the Kish sample, AC2 and for 
HOPG-3 the measured values are larger than the max- 
imum possible contribution from Fe clusters. Thus, for 
these three samples the values of magnetization cannot 
be explained by Fe impurities. 

Moreover, the assumption that such a low Fe- 
concentration behaves ferromagnetically and produces a 
hysteresis loop with finite remanent moment and coercive 
fields at room temperature or even higher 0] is rather 
unlikely taking into account that the Fe concentration 
measured in the HOPG samples was found to be ho- 
mogeneously distributed. To quantify the magnetization 
behavior of Fe impurities in graphite we have measured a 
natural graphite bulk sample that has an inhomogeneous 
Fe distribution through the sample between 13 /Ug/g and 
3,790 ^g/g; most of the sample shows a concentration 
larger than 500 ng/g Fe. All other magnetic impurities 
are below 10 ^tg/g and may therefore be neglected. Fig- 
ure ||(a) shows the hysteresis loop of the natural graphite 
crystal for fields parallel to its main plane after subtract- 
ing the paramagnetic (10 K) or diamagnetic (300 K) con- 
tributions (proportional to the applied field) shown in 
the inset. The measured hysteresis is very small with 
a remanent magnetization M r < 2 x 10~ 5 emu/g. The 
magnetization at 2kOe is < 4 x 10 -3 emu/g. Naively, one 
would argue that such a magnetization can be produced 
by ~ 18 /ig/g ferromagnetic Fe. Such a comparison is, 
however, misleading. If the magnetic signal is due to Fe 
impurities, the nearly negligible hysteresis width and the 
lack of saturation at 2 kOe indicate that Fe impurities 
do not behave ferromagnetically but rather superpara- 
magnetically. To better understand this point we mea- 
sured the magnetization up to 5 T and we obtained a 
ferromagnetic-like curve of negligible width, see Fig. ^(b) , 
after subtracting the paramagnetic background contri- 
bution. We estimate that the saturation magnetization 
observed at 5 T is produced by 840 /^g/g Fe, a concentra- 
tion comparable to the average impurity content of the 
sample. 

A typical response of superparamagnetic Fe in graphite 
is also recognized in the strong temperature dependence 
of the magnetization at 0.2 T below 150 K, see inset in 
Fig. ||(b). This behavior is in clear contrast to that ob- 
tained in the HOPG samples. As an example we show 
the temperature dependence of the magnetization mea- 
sured at 2kOe in the HOPG-2 sample (after subtracting 
the diamagnetic contribution) in Fig. [h]. It is recog- 
nized that the magnetization is practically temperature 
independent, in clear contrast to the expected behavior 
if the 8 fig/g of Fe would behave superparamagnetically. 
This result agrees with that reported in Rcf. Q that the 
saturation magnetization in annealed HOPG samples re- 
mains practically temperature independent up to 800 K. 
In line with that, the measured hysteresis loops for the 
HOPG-2 sample at 800 K (which is the highest possi- 
ble in our equipment) indicate that the Curie temper- 



ature is larger than that. From a linear extrapolation 
of the remanent magnetization (Fig. [l0|) to zero we esti- 
mate that the Curie temperature is of the order of 900 K. 
The small change observed in the hysteresis upon anneal- 
ing the sample at 700 K or 800 K in vacuum, shown in 
Fig. |i|(b), appears to be difficult to explain by arguing 
in terms of a magnetic contribution of ferrocene. From 
all the evidence presented here we believe that magnetic 
impurities can be excluded as the origin of the weak fer- 
romagnetism observed in HOPG . 

(b) Topological defects: Grain boundaries and edge 
states. The topological and electronic structure of atom- 
ically flat graphite surfaces has been extensively studied 
by scanning tunneling microscopy (STM) 
More recently, the electronic properties of carbon nan- 
otubes, and particularly edge states of nanometer-sized 
graphite ribbons, have attracted much interest [li|. It 
is expected that the electronic structure of the graphite tt 
system may differ considerably on carbons located within 
a graphite sheet and those at its edge [^0|, |2l]| . More- 
over, such localized edge states have been shown to de- 
pend sensitively on the geometry of the edge (l8| : While 
"zigzag" edges perpendicular to the [1,0,1,0] lattice di- 
rection show such states, "armchair" edges (perpendicu- 
lar to [1, 1, 2, 0]) do not. Thus, "zigzag" edge states may 
lead to an increase of the density of states at the Fermi 
level. This has been experimentally confirmed in STM 
work on HOPG samples p9| . If such edge states occur at 
high density, a ferrimagnetic spin polarization may result. 
Experimental findings on graphitized nanodiamond [22^ ] 
and activated carbon fibers p3[ support this expectation. 
Similarly, a strong enhancement of the Curie-like param- 
agnetic contribution was predicted at low temperatures, 
see Ref. jlOj and references therein, and this has been ex- 
perimentally confirmed pjj ]. Recently, the mechanisms 
of magnetism in stacked nanographite were theoretically 
studied f25f . This work obtained an antiferromagnetic 
solution for A-B-type stacking only. 

One may doubt that the edge states could be the origin 
for the ferromagnetism found in our macroscopic graphite 
samples, since in nanographites these states provide an 
enhanced paramagnetism and their size ~ 2 nm is smaller 
than the in-plane correlation length of our graphite sam- 
ples (the influence of the edge states to the electronic 
properties of a graphene layer should decrease with D). 
However, small-size effects and the thermal energy could 
also preclude the formation of a stable ferromagnetic or- 
dering in nanographite. It is also unclear whether corre- 
lation effects between edge states in macroscopic graphite 
samples with a in-plane correlation length not more than 
five times the size of the reported nanographite samples 
(see Table I) would not be enough to stabilize a ferromag- 
netic order. Theoretical studies suggest that the pres- 
ence of strong topological disorder will stabilize ferro- 
magnetism and frustate antiferromagnetic order No 
prediction was made, however, on the minimum density 
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of the disorder necessary for this situation . Therefore 
and because the theoretical results are not yet conclu- 
sive we should not neglect a possible contribution of edge 
states to stabilize the FM in bulk graphite. 

In order to assess the role of edge states and their 
potential contribution to the magnetic properties of the 
graphite samples studied in this work, a systematic STM 
investigation of these samples would be desireable. Un- 
fortunately, such a study is beyond the scope of our cur- 
rent experimental capabilities. In order to obtain at least 
a feeling whether or not edge states that might exist at 
grain boundaries within the graphite samples contribute 
to their observed anomalous magnetic properties we have 
conducted a limited study of some of the samples in air. 
We note that we are currently not capable of perform- 
ing the spectroscopic studies [^6| that were required to 
assess quantitatively putative differences in the density 
of states at the Fermi level as a function of the atom 
positions with respect to grain boundaries. 

Figure [ll] gives exemplary STM results obtained at 
a freshly cleaved surface of the HOPG-2 sample in air. 
The images have been flattened, but no filtering was ap- 
plied to the raw data. Panel (a) shows a 125 x 125 nm 2 
surface topology overview observed in constant-current 
mode. Two types of linear defects are observed within 
the planar surface: A step line, oriented roughly along 
the scan direction (horizontal direction within the im- 
age) and an apparent elevation line that runs along the 
step defect at an oblique angle. Both line defects are 
much longer than the width of the image and are char- 
acterized in the section profile shown on the right. In 
other micrographs than the one shown in Fig. [ll] we do 
not find indications that the two type of line defects are 
coupled to each other. By inspection of a large number 
of similar micrographs, we conclude that such defects are 
characteristic of the sample surface. On the surface of 
the HOPG-2 sample, they have a typical distance of a 
few 100 nm. Inspection of the unit cell orientations on 
both sides of the two types of defects at higher resolution 
(vide infra) reveals that both are associated with grain 
boundaries. 

The section profile indicates that the step height of the 
first line defect is ~ 0.65 nm, which corresponds to the 
length of one unit cell along the crystallographic c axis, 
i.e. the height of two graphite sheets. In fact, an addi- 
tional feature along the defect line at mid- level of the step 
is recognized. Since two stacked it rings in the unit cell 
are displaced sideways with respect to each other in the 
graphite crystal, this half-step may correspond to an indi- 
vidual graphite layer that protrudes laterally underneath 
the top layer. At a distance of a few nm from the step 
line, the profile is essentially flat. This is expected as the 
two planes correspond presumably to equivalent crystal 
surfaces that differ by one unit cell in height. However, as 
the trace on the higher level approaches the step, a pro- 
nounced steepening of the slope is observed. Since the 



scan direction was chosen to be parallel to the step line, it 
cannot be argued that this feature is due to tip overshoot. 
Rather, because the experiment is inherently incapable of 
discriminating between topological features and changes 
of the density of states, the observed steepening of the 
contour line may more likely be associated with an in- 
crease of the local density of states on the carbon atoms 
as the topology approaches the grain boundary. 

The apparent topology of the second line defect is dis- 
tinctly different. The section profile registers a sharp, 
symmetric peak that has an apparent height, Az ~ 
0.25 nm, and an apparent width of w 2 nm. At higher 
resolution (Fig. |ll|(b) , showing a constant-height scan of 
a similar defect on the same sample surface) the defect 
structure is disclosed in more detail. As recognized in the 
overview (inset) the defect line appears as three collinear 
stripes. These are no tip artifacts since they do not de- 
pend on the relative orientation between the defect line 
and the scan direction. The high resolution micrograph 
reveals that the corrugation is associated with individual 
surface atoms and that the lattice orientation differs on 
both sides of the central spine. The apparent height of 
the features with the highest elevations - truncated in the 
grey level scale which was chosen to better visualize the 
atomic lattice on the surrounding crystal faces - is about 
a factor of two larger than an atomic diameter. It is thus 
likely that an increase in the density of states on atoms 
near the grain boundary contribute to the appearance 
of these corrugations. Without a clean discrimination 
of topological and density-of-state features, however, it 
cannot in this case be decided with confidence whether 
or not the observed defects correspond simply to lines of 
atoms that decorate the grain boundary. 

A high-resolution characterization of step-type line de- 
fects is presented in the lower half of Fig. O. Panels 
(c) and (d) show the vicinity of the defect shown in 
Fig. [ll] (a) at high resolution in constant-current and 
constant-height mode, respectively. The scan direction 
has been rotated by « 90° with respect to that in the 
overview scan. Similar to Fig. [ll] (b), the images show 
the surface at atomic resolution. The fundamental cor- 
rugation is due to a carbon site asymmetry JTJ, [l5) giv- 
ing rise to a hexagonal lattice. Again, the constant- 
current micrograph - Fig. |ll](c) - indicates a step height 
of Az » 0.6 nm (section profile on the right). Both mi- 
crographs reveal clearly a misorientation of the lattices 
on both sides of the defect line of 42° ± 5°. On the higher 
level - to the left - the main lattice direction meets the 
grain boundary at 65° ± 5°, characteristic of a "zigzag" 
edge. By contrast, one would expect this angle to be 90° 
for an "armchair" edge. In distinction from an earlier re- 
port [|l9) , where "zigzag" edges were observed to extend 
only for short lengths, the defects observed here may ex- 
tend for > 100 nm, as demonstrated in Fig. [ll] (a) where 
the step defect runs along the same general direction over 
the entire, 125 nm, width of the micrograph. 
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In the close vicinity of the grain boundary itself, the 
constant-height image, Fig. |ll] (d), appears distinctly dif- 
ferent from the constant-current image. It is character- 
ized by a vast overshooting of the section profile at the 
locations of the edge carbon atoms, c.f. pseudo-3D rep- 
resentation on the right of Fig. |Tl|(d). While we reiter- 
ate that a clean distinction of topological and density- 
of-state features is not afforded in our measurements, we 
may note that the observed effect cannot easily be ac- 
counted for by the sample topology or elasticity effects 
p6| . For example, we consider it unrealistic that a deco- 
ration of the grain boundary with atoms (of an ordinary 
- low - density of states at the Fermi level) should be 
responsible for the observed giant corrugations. Rather, 
a large increase of the local density of states at the atom 
positions near the edge must be invoked to explain the 
observations. This is also hinted by the deviation of the 
section profile in the overview scan, Fig. [ll] (a), from a 
planar topology near the step. The conjectured enhance- 
ment of the density of states on edge atoms may indicate 
the existence of a localized state as demonstrated in sim- 
ulations [ p8| . Therefore, one may argue that topological 
defects such as the "zigzag" edge characterized in Fig. [Tl| 
affect the electronic properties of the bulk graphite crys- 
tal and might also be responsible, at least partially, for 
the measured ferromagnetic behavior. Clearly, a more 
definite conclusion on the importance of such localized 
states requires much more systematic STM work. In 
future, we are also planning to vary systematically the 
number of topological defects and thus quantify their in- 
fluence on the magnetization. Preliminary results indi- 
cate indeed an increase of the ferromagnetic signal upon 
milling of a bulk sample. 

Finally, we note that recent theoretical work points 
out the effects of the electron-electron interactions in a 
graphene layer suggesting the existence of both ferromag- 
netism and an anisotropic p— wave superconducting state 
||. This work also suggests that topological disorder in 
the graphene planes may enhance the electronic density 
of states and induce instabilities in the electronic system, 
thereby giving rise to the behavior reported here and in 
earlier work 

(c) Itinerant ferromagnetism. Since the carrier den- 
sity in graphite is low, the magnetic properties may 
be connected to itinerant ferromagnetism of dilute two- 
dimensional electron gas systems when the electron- 
electron interaction is large. Indirect evidence for the im- 
portance of the electron-electron interaction in graphite 
has been recently obtained by magnetoresistance mea- 
surements that showed a field induced metal-insulator 
like transition with a scaling quantitatively similar to 
those found in MOSFET's and amorphous thin films [p9| . 

In general, the strength of the electron-electron inter- 
action can be estimated through the Coulomb coupling 



constant 



(Trn 2 D) 1/2 a* l 



(1) 



where the 2D carrier density is given by n2D = n^od with 
the interplane distance d — 0.335 nm, the effective Bohr 
radius a* B — eh 2 /e 2 m*, and the dielectric constant for 
graphite e = 2.8. The 3D carrier density of the majority 
band is n^D — 2 x 10 18 cm -3 with an effective mass 
to* ~ 0.05 mo where too is the free electron mass. For the 
minority carriers, n^u — 6 x 10 16 cm~ 3 with an effective 
mass to* ~ 0.004 to [|o|. Substituting these values in 
eq. ([!]) we obtain r s ~ 5 and ~ 10 for the majority and 
minority carriers, respectively. A value of r s larger than 
1 implies that the potential energy per electron is larger 
than the Fermi energy, hence effects due to the electron- 
electron interaction must be taken into account. 

If we assume that Dirac fermions are responsible for 
the magnetic field driven metal-insulator-like transition 
in graphite (for fields parallel to the c— axis), the strength 
of the Coulomb coupling is characterized by the massless 
and dimensionless parameter pll B2j 
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eQVph ' 



(2) 



where e is the electronic charge, eo the dielectric constant 
and vp ~ 2 x 10 6 m/ s the Fermi velocity for graphite. 
As noted recently [pij , for graphite, g > 10 holds which 
suggests that the strong Coulomb interaction can open 
an excitonic gap in the spectrum of Dirac fermions. In- 
terestingly, the occurrence of the gap is accompanied by 
the appearance of a small magnetic moment as the result 
of band anisotropy |3lf| . In agreement with this theory, 
conduction electron spin resonance (CESR) experiments 
performed in HOPG and Kish graphite samples reveal a 
ferromagnetic-like internal field J|j|. It is expected that 
this doped excitonic FM would disappear above a certain 
level of doping of the order of the chemical potential |3l) . 
We note that within this picture the excitonic gap and 
the nonzero spin polarization is driven by the Coulomb 
interaction even at zero magnetic field. On the other 
hand, a field applied normal to the graphene layers can 
induce the formation of a gap even when the Coulomb 
interaction is weak pi] ], a model used to interpret the 
recently found metal-isolator transition in HOPG p9[ . 

It is also interesting to compare the magnetism found 
in graphite with that in the high-temperature itinerant 
Cai^zLa^Bg ferromagnet |34| where either quantum lo- 
calization effects |35| or an excitonic mechanism [[36] arc 
invoked to explain the observed weak ferromagnetism 
with a Curie temperature of 600 K . In graphite the low- 
temperature ferromagnetic magnetization at saturation 
is ~ 1 x 10 -3 emu/g. This value can be translated in 
terms of the magnetic moment ~ 0.1 hb per majority 
carrier similar to 0.07/zb found in Cai-ajLa^Bg. Accord- 
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ing to Ref. |3lj the ferromagnetism in graphite and hex- 
aborides can have the same origin. Direct experimental 
evidence supporting an excitonic mechanism for the FM 
found in HOPG is, however, still lacking. 

5. Summary 

We have studied the magnetization of several HOPG, 
Kish graphite and natural graphite samples with dif- 
ferent content of magnetic impurities. Our results rule 
out that ferromagnetic impurities can be responsible 
for the observed magnetic effects in HOPG and Kish 
graphite. The magnetization at saturation remains prac- 
tically temperature- independent up to 500 K. This be- 
havior is in clear contrast to that found in natural 
graphite samples with at least two orders of magnitude 
larger concentration of Fe-impurities. The results of nat- 
ural graphite reveal the superparamagnetic behavior of 
Fe impurities in the carbon matrix. Long-time annealing 
at 700 K and 800 K in vacuum of the HOPG samples 
affects the hysteresis loops only weakly. 

The origin of the FM in HOPG is not yet clear. We 
have discussed two possible origins in this work. Topo- 
logical defects can contribute to the peculiarities of the 
graphite electronic properties and may give rise to fer- 
romagnetic correlations. Also, the strong Coulomb in- 
teraction between electrons in graphite (due to the small 
electronic density) should play an important role in all 
the observed magnetic properties of graphite. However, 
details of the long-range itinerant ferromagnetic order re- 
main unclear. 

The discovery of ferromagnetism in oxidized Ceo §J 
as well as in polymerized rhombohedral Cgo indicates 
that carbon-based materials can be ferromagnetic with- 
out metallic components at high temperatures. There- 
fore, it should be no surprise if graphite itself shows weak 
ferromagnetic signals. Future work should try to (a) en- 
hance the ferromagnetic signals by an appropriate mixing 
of graphite with other elements (an enhancement of the 
superconducting signal has been recently achieved mix- 
ing graphite with sulfur jll], |37j ) , (b) increase the defect 
concentration as well as (c) prove experimentally the ex- 
istence of an excitonic gap. 
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TABLE I: Characteristics of the investigated samples. The samples AC1, AC2, UC3, UC4, HOPG-2 and HOPG-3 are highly 
oriented pyrolytic graphite (see text), NG is a natural graphite crystal from Ticonderoga (USA) J^J. (a) All magnetization 
values are in units of 10 -4 emu/g. The magnetization values in the table were obtained after subtraction of a linear background 
contribution, (b) The Fe concentration in this sample is not homogeneous. In four different position of the sample the Fe 
concentration was found to be between 0.10% and 0.38% (wt), in two other positions the values were between 13 /ig/g and 165 
/ig/g. (c) All impurity concentration values are in /ig/g. For all samples the MDL < 2 /ig/g. n.d.: not determined. 
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FIG. 1: (a) Magnetic moment of the sample holder (plastic tube with a small window and an adhesive tape) as a function of 
applied field at T = 10 K. (b) The results from (a) after subtraction of the linear paramagnetic background. 1 emu/g = 1 A 
m 2 / kg, 1 Oe = 10 3 /4?r A/m. 
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FIG. 2: (a) Magnetization as a function of applied field at T — 300 K for sample AC1. No background subtraction has been 
performed, (b) The same for sample AC2 at 10 K after subtraction of a diamagnetic background, \ = —1.32 x 10 -6 emu/g 
Oe. The sample holder background (see Fig. hi) was not subtracted from the data. 
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FIG. 3: (a) Magnetization as a function of applied field at T = 300 K for sample UC3 (Union Carbide). A diamagnetic 
background \ = —2.58 x 10~ 6 emu/g Oe has been subtracted, (b) The same for sample UC4 at 300 K after subtraction of a 
diamagnetic background \ = —1-74 x 10~ 6 emu/g Oe. 
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FIG. 4: (a) Magnetization as a function of applied field at T = 10 K for sample HOPG-2. (b) Similar plot for a different piece 
of the same sample at 300 K after annealing in vacuum at a temperature of 700 K for 16 h (o) and after 19 h at 800 K (▲). The 
bold line represents the hysteresis at the same temperature before annealing. A diamagnetic background has been subtracted. 
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FIG. 5: Remanent magnetization as a function of Fe concentration at two temperatures of all HOPG and Kish graphite samples. 
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Fe Concentration (ng/g) 

FIG. 6: Magnetization at 2 kOe as a function of Fe concentration at two temperatures of HOPG and Kish graphite samples. 
The solid line represents the expected magnetization if Fe contained in the samples was in a ferromagnetic state. The dashed 
line represents a similar relation for FesCU. 
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FIG. 7: Coercive field at two temperatures of the HOPG and Kish graphite samples as a function of Fe concentration. 
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FIG. 8: (a) Magnetization at 2kOc (T = 10 K) as a function of the full width at half maximum (FWHM) in 9 for all HOPG 
and Kish graphite samples, (b) Similar plot as a function of the FWHM in 29. 
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FIG. 9: (a) Magnetization as a function of applied field at 10 K and 300 K for the natural graphite crystal with an average 
Fe-concentration of more than 500 /Ug/g. A linear background contribution was subtracted. These are shown in the inset, (b) 
Similar plot at 10 K and up to ±50 kOe applied field. The inset shows the temperature dependence of the magnetization at 2 
kOe. 
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FIG. 10: Temperature dependence of the saturation (■) and remanence (□) magnetic moment for the sample HOPG-2 after 
annealing 16 h at 700 K. In the present work and due to technical reasons we were not able to measure M(T) continuously up 
to 800 K. 
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FIG. 11: STM images of a freshly cleaved surface of the sample, HOPG-2, in air (left) and their analysis using section profiles 
(right). The ordinate scales in (a) - (c) apply to both the section profiles and the grey scales, (a) Overview scan in constant- 
current mode (Vt = 32 mV, Io = 1.25 nA, / = 6 Hz). In the section profile determined along the dashed line (right), a 
linear ramp has been subtracted to emphasize the flatness of the plateaus. (Thus, the profile levels do not entirely faithfully 
correspond to the grey levels shown in the micrograph.) Close to the step line that is recognized in the lower third of the 
image, an increase of the slope of the section profile is observed. Correcting for this overshoot, a step height of w 0.6 nm 
is determined, (b) High- resolution constant-height image and cross-section (Vt = 33.5 mV, Jo = 1.0 nA, / = 27.5 Hz) of a 
surface region incorporating a line defect similar to the one shown in panel (a) "north" of the step defect. Atomic resolution is 
achieved within the planar sections of the sample surface and indicates that the line defect corresponds to a boundary between 
differently oriented grains. The inset in the micrograph shows a constant-current (Vt = 33.5 mV, Io = 1.0 nA, / = 2.0 Hz) 
overview of the same surface feature. Scan size here is 25.4 x 25.4 nm . The grey scale of the overview is distinct from that 
of the main image and ranges from to 1.0 nm (black to white), (c) and (d) Constant-current (Vt = 20 mV, Io = 2.0 nA, 
/ = 2.9 Hz) and constant-height (Vt = 20 mV, Io = 2.0 nA, / = 55 Hz) images, respectively, of the step defect shown in panel 
(a). The scan direction has been rotated by 90° with respect to that in the overview scan. The section profile in (c) has been 
determined in the rectangle bounded by thin dashed lines. Fat dashed lines indicate the grain boundary as well as the prevalent 
lattice directions within the grains near the boundary. In panel (d), the image is dominated by giant corrugations at atom sites 
immediately at the boundary. The relative proportions of these corrugations with respect to the amplitudes observed further 
away from the defect line is more clearly revealed in the pseudo-3D plot shown on the right. 



